Electronic structures of III-V zinc-blende semiconductors from atomistic first principles 
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For analyzing quantum transport in semiconductor devices, accurate electronic structures are critical for quan- 
titative predictions. Here we report theoretical analysis of electronic structures of all III-V zinc-blende semicon- 
ductor compounds. Our calculations are from density functional theory with the semi-local exchange proposed 
recently [F. Tran and P. Blaha, Phys. Rev. Lett. 102, 226401 (2009)], within the linear muffin tin orbital scheme. 
The calculated band gaps and effective masses are compared to experimental data and good quantitative agree- 
ment is obtained. Using the theoretical scheme presented here, quantum transport in nanostructures of III-V 
compounds can be confidently predicted. 
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I. INTRODUCTION 

During the past five decades, semiconductor device minia- 
turization has brought modern device technology to the 
nanometer scale where quantum phenomena of charge and 
spin transport start to dominate the device physics. 1,2 For tiny 
devices, the atomic nature of the materials are playing an in- 
creasingly prominent rolei and charge transport in these sys- 
tems driven by external fields is intrinsically a nonequilib- 
rium problem. As a result the operation of nano-scale elec- 
tronics crucially depends on the close coupling of nonequilib- 
rium quantum transport phenomena with the atomic structure 
of the device material. Such a situation poses serious chal- 
lenges to theoretical understanding and computational model- 
ing of nanoelectronic device physics. More recently, atomistic 
methods have been combined with the Keldysh nonequilib- 
rium Green's function (NEGF) formalism to meet the chal- 
lenges of quantitative analysis of nanoelectronics. The atom- 
istic methods are used to determine the material properties as 
well as the device Hamiltonian while NEGF is used to pre- 
dict the nonequilibrium density matrix and transport proper- 
ties. When combined self-consistently, such techniques can 
predict not only qualitative but also quantitative properties 
of quantum transport in nano-devices. Along this line, the 
state-of-the-art formalism is to carry out self-consistent den- 
sity functional theory (DFT) atomistic calculation within the 
NEGF framework.- 5 So far, parameter-free NEGF-DFT meth- 
ods have been widely applied to capture quantum transport 
physics from the atomistic point of view4. 

However, in order to apply first principles approach to in- 
vestigate semiconductor nanoelectronics, some very serious 
issues have to be resolved. First, realistic semiconductor de- 
vices (e.g. transistors) have large number of atoms and doped 
with small concentrations of impurities, while typical DFT 
methods can only comfortably deal with low hundreds of 
atoms. Second, DFT with the local density approximation 
(LDA)i2r— and generalized gradient approximation (GGA)^ 



underestimate band gaps of semiconductors. One could not 
predict transport results if band gaps and dispersions were 
not accurate: this is especially serious for semiconductors ap- 
pearing in transistors since their gaps are not very large to 
begin with. Advanced methods such as GW— and hybrid 
functional^ can yield accurate band gaps for many systems, 
but require very large computation for semiconductor devices 
having hundreds to thousands of atoms. 

For pure semiconductors, a recently proposed modified 
Becke-Johnson (MBJ) semilocal exchange was shown to give 
good band gaps for many semiconductors with a computa- 
tional cost similar to that of MBJ is not a fundamental 
solution to the issue of electron correlation, but it is practically 
useful for calculating band structures thus helpful for analyz- 
ing quantum transport properties. By implementing the MBJ 
functional into a newer generation of NEGF-DFT technique 
which is based on linear muffin tin orbital (LMTO) method 
with the atomic sphere approximation (ASA),- transport in 
Si nano-transistors with localized doping and large number of 
atoms was recently analyzedJi Clearly, a very important next 
step is to investigate III-V semiconductors. 

The III-V compound semiconductors are the most im- 
portant materials for optoelectronic device applications and 
are also very important for the complementary metal-oxide- 
semiconductors (CMOS) technology. Major efforts are 
devoted by the microelectronics industry to integrate III- 
V semiconductors into Si CMOS^ III-V semiconductors 
have been extensively investigated both experimentally and 
theoreticallyr^H with particular attention paid to their band 
topologies since band parameters are critical for understand- 
ing quantum transport. Electronic structures of several direct 
gap III-V compounds have also been calculated using the MBJ 
functional as implemented in the planewave DFT codesJS It 
is the purpose of this work to employ the DFT-MBJ approach 
within the LMTO-ASA scheme, to accurately calculate the 
band parameters of all the zinc-blend III-V semiconductors. 
Our calculated band gaps at high-symmetry points (r, X, and 
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L) are quantitatively compared with the corresponding exper- 
imental data; our calculated effective hole/electron masses at 
r point are also well consistent to the experimental values. 
Details of the ASA schemes and the MBJ potentials will also 
be presented. 

The rest of the paper is organized as follows. In the next 
section the calculation method is briefly discussed. Section 
III presents the results and Section IV is a short summary. 



II. METHOD 



TABLE I: The lattice constant a of the III-V compounds and the c- 
values used in the DFT-MBJ calculations. c a t om and c vac are the 
c-values for real atoms and the vacancy spheres, respectively. 



Material 


a(A) 


Catom 


Cvac 


Material 


o(A) 


Catom 


Cvac 


GaAs 


5.6533 


1.20 


1.39 


GaP 


5.4505 


1.13 


1.50 


AlAs 


5.6611 


l.ll 


1.44 


A1P 


5.4672 


0.69 


1.44 


InAs 


6.0583 


1.18 


1.00 


InP 


5.8697 


1.04 


1.39 


GaSb 


6.0959 


1.17 


1.23 


GaN 


4.50 


1.19 


1.50 


AlSb 


6.1355 


1.12 


1.33 


A1N 


4.38 


1.59 


1.55 


InSb 


6.4794 


1.19 


0.62 


InN 


4.98 


1.56 


1.39 



Our calculation is based on the DFT-MBJ self-consistent 
approach where DFT is within the LMTO scheme and the 
atomic sphere approximation^ as implemented in the Nan- 
odsim software packaged For technical details of the Nan- 
odsim algorithm we refer interested readers to the original 
literatureJ^l A 12 x 12 x 12 fc-mesh were used to sample 
the Brillouin Zone of the primitive cell. The lattice constant 
for the semiconductors were adopted from Ref. [l6| and listed 
in Table U In order to carry out LMTO DFT calculations for 
semiconductors, a good ASA scheme is very helpful. In our 
ASA, vacancy spheres were placed at appropriate locations as 
done in Ref. [22j for space filling, and the same sphere radius 
are used for all the vacancy spheres and atomic spheres. Elec- 
trons in the full d-orbital of Ga, In, As, and Sb are included as 
valence electrons. After the LMTO DFT self-consistent calcu- 
lation is completed, band structures are calculated by the the 
muffin-tin orbital (MTO) approach. The effective hole mass 
m* h and effective electron mass to* were then obtained by fit- 
ting the valence band maximum and conduction band mini- 
mum to a parabola, respectively. Spin-orbit coupling was not 
considered in this work. 

Some details of the MBJ semilocal exchange potential is 
worth mentioning. The MBJ potential has the following 
form^ 2 - 



«& BJ (r) = c <«(r) + (3c-2)ijA .ft^, (1) 

7T V 12 y p a (r) 

where subscript a is spin index, p a is the electron density, 
t a is the kinetic energy density, and v^(r) is the Becke- 
Roussel potential, 23 The relative weight of the two terms is 
given by a parameter c which depends linearly on the square 
root of \\7p\/p. For all solids investigated by the MBJ poten- 
tial so far^ it appears that E g increases monotonically with 
c. The c-value can be determined by the protocol discussed 
in Refs.dUl In our calculations, because the LMTO- ASA 
is a site oriented technique, it allows the c-value to be "local", 
namely one can use a c-value for the atom and another c-value 
for the vacancy sphere. In particular, for a given compound, 
the same c-value is used for the real atoms (e.g. Ga and As), 
and another c-value is used for the vacancy spheres. Recently, 
this MBJ scheme was applied in Ref.|25]to determine very ac- 
curate band gaps for the Al^Gai-^As compounds and band 
offsets for the GaAs/Al 2 .Ga 1 _ 2 .As heteroj unctions for the en- 
tire range of the concentration < x < 1. Here we use the 
same scheme for all the III-V compounds and the optimized 



c-values are found and listed in Table J] Even though one may 
expect more accurate band parameters by using different opti- 
mized c-values for all the different atomic spheres, the scheme 
we use is a good compromise between being simple and also 
reasonably accurate. 

III. RESULTS 

We begin by calculating the band structures of the III-V 
compounds with LDA 5 using the projector augmented waves 
(PAW) method and a plane-wave basis set of 400 eV, as imple- 
mented in the Vienna ab initio simulation package (VASP); 26 
as well as using the LMTO-ASA method discussed in the 
last Section where ASA is given in Ref. |22j, as implemented 
in the Nanodsim 13 NEGF-DFT package. Results for all the 
zinc -blende III-V compounds are plotted in Fig. (fTJ, show- 
ing a perfect agreement of the valence bands and a very good 
agreement of the conduction bands between these methods, 
confirming that our ASA scheme in the LMTO-ASA is accu- 
rate for calculating physical properties of these materials. It 
should be noted from Fig. (fl3 that the band gaps were under- 
estimated by LDA. 

Having confirmed our ASA scheme, we next apply the MBJ 
functiona l 12 ' 24 to calculate the electronic structure again using 
the LMTO-ASA approach, and the MBJ results are plotted 
in Fig. ©. Compared with the LDA bands in Fig. ([TJ, the 
opening of band gaps is evident. From Fig. (O, the band gap 
values at T, X and L points are obtained and listed in the 
third column of TableHII In Ref.yjl five direct gap compounds 
were calculated by the MBJ functional within a planewave 
method, their band gap values are also listed in Table Alias the 
fourth column. Our results for these five compounds agree 
very well with the planewave results except for the X point 
of InAs and InSb. As discussed in Ref. [Tgj, for these narrow 
band semiconductors, it is difficult to accurately determine the 
gaps at the X point. For instance, exp erimental gap value of 
0.63 eV was recommended in Ref. [Ta but other values such as 
1.80 eV was also reported in earlier literature^. 

The experimental gap values are taken from Ref. Il6lll7l and 
listed as the fifth column in Table [XT] The last column is the 
percentage difference between our calculated values and the 
experimental values. Of the 36 tabulated gap values, 32 of 
them are in very good agreement with the experimental data. 
There are 4 values having a difference over 20%: the L point 
of InSb; the X and L points of InN; and the L point of AlSb. 



3 



TABLE II: Energies of the conduction band minima (band gaps, E g ) 
at the r, X, and L points with respect to the valence band maxi- 
mum at the r point in units of eV, calculated by the DFT-MBJ ap- 
proach at zero temperature. The experimental results are taken from 
Ref. [If] (except where noted). MB Jv shows the DFT-MBJ values 
from Ref. El The minus sign of the deviation indicates that the cal- 
culated value is smaller then the experimental value. 



Material 




This work MBJy " 


Expl. 1 " 


Deviation(%) 


GaAs 


r 


1.529 


1.52 


1.519 


0.7 




X 


2.003 


2.00 


1.981 


1.1 




L 


1.682 


1.72 


1.815 


-7.3 


AlAs 


r 


3.087 




3.099 


-0.4 




X 


2.240 




2.240 







L 


2.800 




2.460 


13.8 


InAs 


r 


0.416 


0.43 


0.417 


-0.2 




X 


1.440 


2.01 


1.433 


0.5 




L 


1.225 


1.43 


1.133 


8.1 


GaP 


r 


2.887 




2.886 







X 


2.350 




2.350 







L 


2.429 




2.720 


-10.7 


A1P 


r 


3.635 




3.630 


0.1 




X 


2.513 




2.520 


-0.3 




L 


3.030 




3.570 


-15.1 


InP 


r 


1.421 


1.42 


1.424 


-0.2 




X 


2.502 


2.34 


2.384 


4.9 




L 


2.007 


2.11 


2.014 


-0.3 


GaSb 


r 


0.818 


0.82 


0.812 


0.7 




X 


1.117 


1.21 


1.141 


-2.1 




L 


0.876 


0.87 


0.875 


0.1 


AlSb 


r 


2.346 




2.386 


-1.7 




X 


1.698 




1.696 


0.1 




L 


1.845 




2.329 


-20.8 


InSb 


r 


0.238 


0.25 


0.235 


1.3 




X 


0.613 


1.52 


0.630 


-2.7 




L 


0.477 


0.82 


0.930 


-48.7 


GaN 


r 


3.298 




3.299 







X 


4.528 




4.520 


0.2 




L 


5.997 




5.590 


7.3 


A1N 


r 


5.853 




6.0 


-2.5 




X 


4.908 




4.9 


0.2 




L 


9.304 




9.3 





InN 


r 


0.781 




0.78^ 


0.1 




X 


3.456 




2.51 


37.7 




L 


4.635 




5.82 


-20.4 



These differences could be due to the uncertainties of the ex- 
perimental values for the narrow gap compounds InSb and 
InNi2 and, of course, the approximative nature of the DFT 
calculations. Nevertheless, it is impressive that, all in all, the 
MBJ band gaps are in very good consistency to the experi- 
mental values for these wide range of materials. 

From the calculated band structures shown in Fig. (|2), we 
obtain the effective masses of electrons and holes by fitting 
the valance band maximum or conduction band minimum to 
a parabola at the T point along the [100] direction, the re- 
sults are listed in Table [HI] together with the experimental 
values-^ 7 The deviation of our calculated values to the ex- 
perimental values are also listed for each compound. Very 
good agreement is obtained for most situations except for A1N 
and the electron effective mass of GaN. Given the fact that ef- 
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FIG. 1: (color online) The band structures for the III-V compound 
semiconductors obtained with LDA. Red line is obtained by VASP, 
blue circles obtained by Nanodsim. 
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FIG. 2: (color online) The band structures for the III-V compound 
semiconductors obtained with MBJ. Blue crosses indicate the exper- 
imental values. 
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TABLE III: Effective hole and electron masses at the T point in units 
of the electron rest mass m e calculated along the [100] direction. The 
experimental values are in the parentheses and taken from Ref. [If] 
(except where noted). The minus sign of the deviation indicates that 
the calculated value is smaller than the experimental value. 



Material 


\ m h\ 


Deviation(%) 


\fn%\ 


Deviation(%) 


GaAs 


0.355(0.350) 


1.4 


0.076(0.067) 


13.4 


AlAs 


0.510(0.472) 


8.1 


0.149(0.150) 


-0.7 


InAs 


0.373(0.333) 


12.0 


0.029(0.026) 


11.5 


GaSb 


0.265(0.250) 


6.0 


0.041(0.039) 


5.1 


AlSb 


0.367(0.357) 


2.8 


0.118(0.140) 


-15.7 


InSb 


0.263(0.263) 





0.018(0.014) 


28.6 


GaP 


0.411(0.326) 


26.1 


0.163(0.130) 


25.4 


A1P 


0.595(0.518) 


14.9 


0.209(0.220) 


-5.0 


InP 


0.450(0.532) 


-15.4 


0.094(0.080) 


17.5 


GaN 


0.927(0.855) 


8.4 


0.232(0.15) 


54.7 


A1N 


1.587(1.020) 


55.6 


0.319(0.25) 


27.6 


InN 


0.977(0.833) 


17.3 


0.085(0.07ii) 


21.4 



fective masses are usually very hard to determine accurately, 
the good consistency to the experimental results for most sit- 
uations is very satisfactory. 

IV. SUMMARY 

We have calculated the electronic structures of all zinc- 
blende III-V semiconductor compounds from density func- 



tional theory with the semi-local exchange of the MB J form, 
using the LMTO-ASA scheme. In our method and due to 
ASA, vacancy spheres are added to fill the volume of the semi- 
conductors. Since this is a sited oriented calculation method, 
the weight parameter c in the MB J potential for the atoms and 
for the vacancy sites are different. We determine the opti- 
mal values of this weight parameter for all the compounds. 
The calculated band gaps are mostly in very good agreement 
with the corresponding experimental data. The obtained ef- 
fective masses are also largely in good agreement to the mea- 
surements. For analyzing quantum transport in semiconductor 
nanoelectronics, accurate electronic structures are very im- 
portant for quantitative predictions. In addition, as we have 
shown recently, the band-offset of heterojunctoins can also 
be accurately predicted using the same method as reported 
here^ Together with the results of this paper, quantum trans- 
port properties of III-V systems can now be confidently cal- 
culated from the atomic point of view since the LMTO-ASA 
method can self-consistently calculate very large number of 
atoms 
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